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During the period of this report research was completed on the A^ 
excited state and low-lying dissociative states of the methoxy radical. 
A manuscript was prepared that reports the characterization of the ^A^ 
electronic state, the excitation energies and Franck-Condon factors for 


the ^A 


_ 2 


E system, and the energies of .intersection between the 


state and the nearby dissociative states. The minimum excitation 
energy needed for predissociation of methoxy is predicted along with 
the corresponding implications for atmospheric chemistry. This 
manuscript was submitted for publication to the Journal of Chemical 
Physics on August 1, 1984. A copy of this manuscript, titled 
“A potential-energy surface study of the ^A^ and low-lying dissociative 
states of the methoxy radical", is attached to this report. It should 
be considered a complete and final report on the phase of the present 
research agreement that pertains to the methoxy radical. 


Travel during this period of the agreement consisted of a trip by 
Dr. Jackels to NASA Langley Research Center in May, 1984. This trip 
was for the purpose of carrying on scientific discussions with Dr. D. 
H. Phillips of NASA pertaining to the funded research. 
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ABSTRACT 


Low-lying bound and dissociative states of the methoxy radical have 
been studied at Cg v conformations using accurate ab initio quantum 
cheaical techniques. A double-zeta quality basis set augmented 
with polarization and diffuse functions has been used throughout. 
The A A^ excited state is found to have an equilibriua conformation 
with Rqq= 1.614 A, Rg^-l.OSl A* and 8 ^ 0 q= 1 O 2 . 6 °, and the three a^ 
harmonic fequencies have been calculated (3239, 1308, and 733 
cm -1 ). The following excitation energies have been estimated: 
l^Aj f X^E vertical excitation, 4.21 eV; IC^A^ 4 X^E vertical 
deexcitation, 3.52 eV; and the *A - X system origin, 3.96 eV. The 
surface intersections of the state with three low lying 

dissociative states ( 4 A 2 » ^ 2 » and 4 E> have been characterized. 

The threshold for predissociation has been estimated as 4.47 eV, 
corresponding to the ^A^ - 4 A 2 intersection. The rate of methoxy 
photolysis is estimated to be too small to warrant inclusion in 
atmospheric models. 

PACS numbers: 31.20.0i, 33.80.Gj, 31.50,+w, and 94. 10. Fa 
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I. INTRODUCTION 

Alkoxy radicals form an important class of hydrocarbon oxidation 
intermediates in both combustion and atmospheric chemistry. Methane, 
which gives rise to the prototypic methoxy radical CH3O, is introduced 
into the atmosphere in very large quantities by natural sources and is 
a global constituent of the troposphere. Although most of the methane 
is oxidized in the troposphere, enough of it reaches higher altitudes 
to give rise to stratospheric methane chemistry. Besides being 
important in the methane oxidation chemistry of both the upper^ and 
lower^ atmospheres, methoxy serves in a limited way as a model for more 
complicated alkoxy radicals. Being a highly reactive species, the CH3O 
radical presents considerable experimental difficulties, but due to its 
relatively small size, it is well suited to accurate theoretical 
investigation. This paper reports the results of a quantum chemical 
study of the methoxy radical, with particular emphasis on its excited 
electronic states. A previous publication' 3 , hereafter referred to as 
I, reported a potential surface study that provided an overview of 
several electronic states of CH3O and a detailed view of the 2 E ground 
state. The present paper reports results that place particular 

4 

emphasis on the excited electronic states of the radical. 

Methoxy is thought to be generated in the atmosphere by the 
reaction of the methylperoxy radical CH3O2 with reduced species such as 
NO, SO2, or HO. Methylperoxy is formed by the reaction of CH3 with 02» 
the methyl radical being the product of hydrogen abstraction from 
methane by either the hydroxyl radical or oxygen atom. This chemistry 
is summarized by the following reactions: 

CH 4 ♦ OH -* CH 3 + H 2 0 ( 1 ) 

n 

ch 3 + o 2 ch 3 o 2 


( 2 ) 
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CH 3 0 2 ♦ HO ■) N0 2 ♦ CH 3 0 (3) 

The major removal mechanism for atmospheric methoxy is thought to be 
its reaction with molecular oxygen to form H0 2 and formaldehyde, 

CH 3 0 ♦ 0 2 -> H0 2 + CH 2 0. (4) 

If the photolysis of CH 3 0 

CH 3 0 » CH 3 ♦ 0 (5) 

were also to occur at a significant rate, it could be an important 
process because it would regenerate the methyl radical, with the 
subsequent formation of CH 3 0 2 and further NO or S0 2 oxidation. This 
cyclical process would have the net effect of destroying 0 2 and 
creating “odd oxygen" species such as N0 2 and 0, which in principle 
could perturb the atmospheric ozone concentration. Consideration of 
Fig. 1, a reproduction of semiquantitative potential curves from paper 
I, suggests that the photolysis process of lowest energy would be 
excitaton to those vibrational levels of the electronic state that 

would be predissociated by the intersecting repulsive states. 

The methoxy radical has been studied by electronic absorption 4 and 

s — 7 ft — 1 n 

emission'' spectroscopy, laser induced flourescence A , laser 

magnetic resonance* , ESR* , and electron photodetachment 
13 14 

spectroscopy. * In recent years, methoxy has been the subject of 

several theoretical investigations using various ab initio 
3 1 5-1Q 

techniques. ' These studies have included characterization of the 

2 E ground and 2 A^ excited states, and have provided barrier estimates 
for isomerization and dissociation. The Jahn-Teller distortion of the 
ground state has been found^' 19 to b e with an energy 

decrease of only about 0.6 kcal/mole. The CH 3 0 spin-orbit coupling has 
also received theoretical attention. 17c, 20 

In this paper are presented the results of a potential energy 



Page 5 


surface study of the X 2 A^ excited state of the raethoxy radical. The 
equilibrium geometry as well as a limited vibrational analysis of this 
excited state species is reported. Transition energies and 
Franck-Condon factors are predicted for both absorption and emission in 
the A^ - ^E system, which is responsible for the main features of the 
CH 3 O visible-UV spectrum. Also, the low-lying repulsive states of 
methoxy, which were given a semiquantitative treatment in I, are 
accurately described in the regions of their intersections with the A^ 
surface. These computed intersection energies provide an estimate of 
the minimum excitation energy required for predissociation. Using this 
estimate, the importance of Rxn. (5) in the atmosphere is evaluated. 

During the course of the present study it was discovered that the 
computer program used to carry out the configuration interaction 
calculations of the 2 E ground state reported in paper I incorrectly 
omitted some double excitations from the master configuration list. 

This error was systematic in that it omitted the same configurations 
for each 2 E conformation but treated the 2 A^ state correctly. The 
effect of correcting this program error, when comparison is made with 
the results of paper I, is a uniform lowering of the extrapolated *E Cl 
energies by approximately 0.010 hartree or 0.27 eV. This has a 
relatively small effect on the shape of the E potential surface, but 
shifts that entire surface relative to the 2 A^ excited state, resulting 
in larger transition energies. The general conclusions of paper I are 
unaffected, but some of the numerical values reported have been 
revised. Accordingly, data tables in the present paper include both 
A-^ and corrected E results. Since the changes are not of a 
qualitative nature, the discussions and analysis of paper I remain 
unchanged and are not reconsidered here. 
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II. GENERAL COMPUTATIONAL METHOD 

The basis set and self-consistent field (SCF) method are the same 
as described in paper I. All computations were carried out at C 3v 
molecular conformations using a Cartesian Gaussian basis set of 
polarized double-zeta quality . 21-22 Equivalence restrictions were 
imposed on the e x and e y molecular orbitals, which constitutes a 
variational constraint for the ^E ground state but not for the ^Aj^ 
state. Relaxation of this constraint allows the ^E SCF energy to 
decrease by about 1 x 10 ^ hartrees; however, at the extrapolated Cl 
level (see below) the decrease is less than 10~ 4 hartrees. 

A different set of computer programs^' 3 was used for the present 
configuration interaction <CI) studies than for those reported in paper 
I. However, the methods are essentially the same. The details of the 
Cl calculations will be discussed in the appropriate sections that 
follow, with only a general description given here. The A k 
perturbation theory technique of Gershgorn and Shavitt^ 4 was used to 
select the most important configurations for inclusion in the wave 
functions. The selection procedure accepted or rejected a particular 
space configuration, with its entire set of spin couplings, based upon 
the sum of the A^ contributions for that set. The configurations were 
first ordered by their estimated energy contributions. Then, beginning 
with the least important one, configurations were rejected in sequence 
until the sum of the energy contributions for the rejected 
configurations was equal to z, the cumulative selection threshold. 

This method is designed to neglect a constant amount of correlation 
energy, and therefore is expected to yield balanced treatments of the 
various electronic states and the different regions of their potential 


energy surfaces. 
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There were two types o£ master configuration lists used in this 
work. In the type of calculation known as "all singles and doubles", 
the selection procedure was based upon a master list containing all 
single and double excitations from the root set configurations. In 
some of the potential surface calculations a "first-order" master 
list^, containing those single and double excitations (with respect to 
the root set) with no more than one electron outside the valence space, 
was generated and subjected to the cumulative threshold selection 
procedure. This treatment is expected to yield balanced results when 
applied to the same state at similar conformations, thus generating 
potential surfaces which are displaced from, but nearly parallel to, 
more accurately determined ones. It is not expected that the 
"first-order" approach will result in such a well balanced treatment 
for different electronic states. Consequently, in this study energy 
differences between states were computed using the "all singles and 
doubles" configuration list. 

To obtain improved Cl convergence and to define meaningful valence 
space orbitals for the first-order calculations, approximate natural 
orbitals (ANO's) were used throughout. At each conformation the ANO's 
were obtained from preliminary Cl calculations using several hundred of 
the most important configurations. After freezing* 4- ' those orbitals 
that were occupied in the root set configurations, the ANO's were 
obtained as the natural orbitals* 6, of these small Cl expansions. 
Throughout this work, the oxygen and carbon Is core orbitals remained 
fully occupied, and the corresponding core correlation ANO's (zero 
occupancy) were discarded. 
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III. FIRST-ORDER POTENTIAL ENERGY SURFACES 

A. Method. 

At each conformation the SCF wave functions and ANO' s were 
obtained for both states by the procedures described in Section II 
above and in paper I. The Cl root sets consisted of the dominant 
configuration for each state 

. . . <3a 1 ) 2 <4a 1 ) 2 <5a 1 ) 2 <le) 4 <2e) 3 < 2 E) (6) 

. . . C3a 1 ) 2 (4a 1 ) 2 <5a 1 ) 1 <le) 4 <2e) 4 ( 2 A 1 ) (7) 

Although in paper I a multiconfiguration root set was used for the 2 E 
state, only the dominant term is needed for the present range of 
important molecular conformations. The and. E first-order master 

lists contained, respectively, 12 159 and 8800 spin couplings (CSF's) . 
Using a cumulative threshold of x=2.5 x 10 * hartree, final 
configuration lists contained approximately 4500C 2 A 1 > and 3400( 2 E> 

CSF's at each conformation. Typically, these truncated configuration 
lists accounted for 98-99* of the estimated (A^> first-order 
correlation energy, and the most important rejected configurations had 
estimated contributions of approximately 7 x 10 - ^ hartree. 

B. Results and discussion 

9 9 

1. ^E and A^ equilibrium geometries 

The C^ v conformations corresponding to the minimum first-order Cl 

9 9 

energies are given in Table I for the ^E and ^A^ states of methoxy. 

For comparison, three other recent ab initio results for the ground 
state are given, and the agreement is reasonable. There is one 
previous theoretical result for the ‘Aj state, which utilized a 
smaller basis set <DZ> in an SCF study. Taking into account the 
differences in method, these results are also in reasonable agreement. 
In paper I a more detailed comparison was made of the same studies for 
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the methoxy ground state. The f transition can be viewed as the 

excitation of an electron from the C-0 tr-bonding orbital (5a^> into the 
nonbonding ir-orbital on oxygen <2e>. In the excited state, the C-0 
bond is formally a one electron bond and the bond length is much longer 
than in the ground state. The C-H bond lengths are shortened a small 
amount by the excitation, and the HCO bond angle is decreased by 6-7°. 
This is consistent with a picture of the -CH 3 moiety gradually assuming 

OQ 

the planar (Dg^) geometry of the free methyl radical CR CH =1.075 A) as 
the C — 0 distance is increased. 

2. Harmonic vibrational energies 

The first-order potential surfaces were least-squares fitted to 
polynomials in the variables 


AR co = 

R C0 ■ 

rS co» 


aR CH = 

r ch * 

r 6 CH' 

(8) 

aS HCO = 

®HCO ■ 

® e HCO' 



where the superscript **e" denotes the value of the coordinate at the 
minimum of the potential function. The general form of the polynomial 
is given by 

E = E 0 * k 200 <AR C0 )2 + k 110 AR C0 AR CH * k 020 <AR CH )2 + k 101 aR C0 a 8 HC0 

+ k 011 AR CH A 9 HC0 + k 002 <A 9 HCO >2 + <9) 

where bond lengths are in bohra, angles in degrees, and the energy in 

hartrees. The full ^E and 2/^ data sets contained, respectively, 86 
and 99 points, which were distributed nonuniformly around the 
equilibrium geometry within the ranges ’ Re C 0±°*2 a o' Re CH i0,1 a o' 

and 0 S HCO - 13°. Quadratic and cubic fits were made to these data sets 
as well as to subsets in which the points furthest from the equilibrium 

geometries had been deleted. The smaller sets contained 59 -t E and 76 

2 

A^ data points. Although the various data fits were fairly similar. 
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the force constants obtained from the cubic fits to the smaller data 
sets yielded the most accurate vibrational frequencies where comparison 
wi-th experiment could be made. The quadratic force constants obtained 
from these cubic functions are presented in Table II. 

The values of the force constant ^200 ^ or ^he two e l ectr °nic 
states are consistent with the C-0 bond being much weaker in the ■ fc A 1 
state than in the ground state. There is considerable variation (20 *> 
in the value of this force constant as obtained from the various 
analytical fits. Since k 2 oo ^as a d° minan t effect on the calculated 
C-0 stretching frequency, attempts were made to reduce the uncertainty 
in this force constant by using more flexible potential functions. A 
full quartic fit to the 99-point ^A^ data set was attempted, but the 
statistical uncertainties in the resulting force constants were 
unacceptably large. To see if a partial quartic fit would be an 
improvement, the cubic polynomial was supplemented by the fourth order 
term in aR^q. Inclusion of this term had little effect on 1<200' an<1 
the fit was unsatisfactory because of numerical difficulties in 
locating the minimum of the potential function. It was concluded that 
the ^A^ data 3et was not complete enough to yield meaningful quartic 
constants. The difficulties involved in selecting an adequate grid of 

points for the purpose of obtaining well-defined higher force constants 

30 

have been discussed by others. 

The Wilson FG formalism was used with the force constants in 

Table II to obtain the harmonic vibrational frequencies and normal 
modes of a-^ symmetry for the CHgO and CD 3 O radicals. These frequencies 
are presented in Table III. For both electronic states, the agreement 
with experimental results is reasonably good. The theoretical results 
of Saebo et al.^ were obtained at the Hartree-Fock level and. 
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according to those authors, should be reduced by lOX before comparison 
with experiment. After this scaling, these theoretical results are 
also consistent with the present ones. Analysis'^ of the 
normal -coordinates in terms of .the internal coordinates gives the 
following dominant contributions: - symmetric C-H stretch; ^2 “ 

symmetric HCO bend (umbrella mode); and ^3 - CO stretch. For most of 
the entries in Table III the internal coordinates are almost entirely 
unmixed, and the above identifications represent nearly all of the 
normal mode character. However, in the ^A^ state of both isotopic 
species the ^3 mode has 18 - 20 JS bending character, with the HCO angle 
decreasing slightly as R^ 0 increases. Also, in the state of CD3O 
the CO stretch and DC O bend are mixed strongly, giving 30 * CO stretch 
character to ^2 and 40 ss b ®nd character to O3. This strong mixing in 
CD3O and its implications for spectral assignments were discussed in 
paper I . 

3 . Anharmonic vibrational energies. 

For comparison with the methoxy spectra, which involve long 
progressions in the C -0 stretch, it is necessary to include anharmonic 
contributions in the theoretical estimates of the higher vibrational 
levels. Two one-dimensional vibrational models were used, both of 
which treated methoxy as a pseudo diatomic molecule. In the first 
approach, the HCO bond angles and H-C bond lengths were frozen at their 
equilibrium values, and the C -0 distance was chosen as the vibrational 
coordinate. The potential curve for this model was defined as the cut 
through the methoxy surface that is generated by varying the (CH3>-0 
distance with the other coordinates fixed. The second approach was 
similar, except that the cut through the potential surface was taken 
along the direction of the ^3 normal coordinate. This coordinate. 
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which in general combines motion along all three internal coordinates, 
is mass weighted, and the pseudo diatomic is treated as a one 
dimensional oscillator with a reduced mass o£ unity. The finite 
difference method described by Tobin and Hinze was used to 
numerically integrate the resulting one-dimensional vibrational 
Schroedinger equations. 

The potential curves for both models were generated from the cubic 
fits to the larger data sets (86 and 99 points) . To get a well behaved 

surface at intermediate values of R^q, it was found necessary to extend 

o 

the ‘‘■A^ data set by including points with aR^q values as large as 
1.0 a Q . Since the cubic potential function, Eq.(9), is designed to 
describe the potential surface near the equilibrium geometry, the cuts 
through this surface were smoothly extended with Morse-type functions 
to larger and smaller values of the vibrational coordinates. In Fig. 2 
are presented the two resulting potential curves for the 2 A^ state. 

They are plotted on the same scale by using mass weighted coordinates 
that in each case would correspond to an oscillator of unit mass. For 
equal values of the mass weighted coordinates, the R^q values for the 
two curves differ slightly as indicated in the figure caption. 

The computed vibrational energy levels, presented in Table IV, 
show that the two anharmonlc models differ very little for the 2 E 
state, because in this case the O 3 normal coordinate is almost pure C-0 
stretch. However, the difference between the two sets of levels is 
much larger for the A^ state where the ^3 normal coordinate includes 
considerable HCO bending character. The normal coordinate potential 
curve is seen in Fig. 2 to rise somewhat less steeply with increasing 
Rcq than does the "C-0 distance" curve, resulting in more closely 
spaced vibrational levels. In the vicinity of the sixth or seventh 
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vibrational level the differences between the two sets of results are 
nearly as large as the intervals between levels. From the electronic 
absorption spectrum (Table I of Ref. 4), the separation of the v '=6 and 

O _ -1 

v'=0 levels of the state is calculated to be 3847 cm . As shown 

in Table IV of the present work, the same separation is calculated to 
be 4286 cm - * using the C-0 distance as vibrational coordinate and 3904 
cm - * using the O 3 normal coordinate. This comparison favors use of the 
normal coordinate to describe the lower lying vibrational levels. A 
similar comparison for the state also favors the normal coordinate 
approach but is less compelling, because in that case the two models 
are not very different. 
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IV. ACCURATE TRANSITION ENERGIES FOR THE A 2 A X - 3f 2 E SYSTEM 
A . Method . 

The procedures used in this part of the study were designed to 
provide well-balanced treatments of the methoxy A and X electronic 
states so that accurate transition energies could be calculated. SCF 
calculations were carried out for each state CEqs. (6)- <7)3 and the 

O 

corresponding ANO's were obtained. Master lists containings 41 740 E 

and 41 047 ^A^ spin conf igurtion3 were generated by including all 

single and double excitations with respect to the appropriate root 

<SCF> configuration. All single excitations were retained, and the 

double excitations were subjected to the A^ selection procedure using a 

cumulative threshold x. Cl energies E(x> were obtained at three values 

of x, fit to a straight line, and extrapolated 34 to zero threshold. 

The smallest threshold used in each case corresponded to inclusion of 

approximately 10 OOO spin configurations. 

To account in an approximate manner for the most important effects 

35 

of quadruple excitations, Davidson's correction formula was used, 

AE q = <1 - C 0 2 )bE d , (10) 

where bE d is the correlation energy due to the entire set of double 
excitations, bE q is the analogous quantitiy due to quadruple 
excitations, and Cq is the coefficient of in the "all doubles” Cl. 

In this work, the quantity bE d was estimated to be equal to the 
extrapolated "singles plus doubles" correlation energy, and Cq was 
taken from the Cl calculation at the smallest of the selection 
thresholds . 

The transition energies have been corrected for the vibrational 
zero point energy. All of the degenerate (e-type> vibrational 

IQ O 

frequencies have been recently calculated for the E state, but a 
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complete A^ set is not available. Therefore, we have neglected the 
zero point contributions of the e-type modes, which is equivalent to 
assuming that they are the same in both states and cancel. For the a^ 
modes the vibrational corrections are explicitly obtained from the 
values reported in Table III. 

B. Results and discussion. 

1 . A^A^ <r X^E vertical excitation energy. 

The vertical excitation energy is taken as the difference between 

the A and X state energies as computed at the ground state equilibrium 

geometry (Table I). For both states, the cumulative selection 

thresholds used in the extrapolated Cl calculations were 40, 25, and 

15 x 10 hartree. In Table V are listed the total energies obtained 

at the various stages of the calculation, and Table VI presents the 

transition energies obtained at these different levels. Inclusion of 

the correlation energy via Cl makes a significant contribution to each 

of the three transition energies. This is to be expected, because the 

C-0 bonding is qualitatively different in the two states. The vertical 

excitation energy is affected very little by either the extrapolation 

procedure or the correction for quadruple excitations. The vibrational 

correction is small but significant and is calculated by including the 

o o 

zero point energies of v ^ _ 3 for the E state and ^-2 * or A i state 

(see Table III). This assumes that only the one normal mode () 3 ) is 
excited in absorption, which is consistent with the experimental 
spectrum. The "best" vertical excitation energy from this work, 4.21 
eV, agrees very well with the position of the intensity maximum of the 
absorption spectrum 4 *, which appears to be one of the vibrational bands 
assigned as v'=4 or 5 (4.23 or 4.31 eV) . This value also agrees very 
well with the position of the most intense band of the fluorescence 
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excitation spectrum, 4.23 eV. 10b 

As a further comparison with experiment, Franck-Condon factors 
were calculated for the A^A^(v') f X^E(v"=0> progression using the 
enharmonic vibrational model with the C-0 distance as the vibrational 
coordinate. It would be preferable to use the vibrational treatments 
based upon the -Og normal coordinates. However, because these 
coordinates are not the same for the two electronic states, a simple 
one dimensional integration of the vibrational overlap is not 
appropriate. Using the coordinate R^q, the Franck-Condon factors as 
relative per centages for v'=0-3 are: 10, 38, 74, 97, 100, 84, 61, 39, 
22, and 11. The largest of these values is at v'=4, which agrees 
reasonably well with the assigned band of maximum intensity (Fig. 1 of 
Ref. 4). However, the Franck-Condon factors calculated for v' greater 
than 4 appear to be somewhat smaller than those reported in Ref. 4. It 
is difficult to make a detailed comparison with the experimental 
intensities because the profile of the spectrum seems to suggest the 
presence of some additional underlying structure at the shorter wave 
lengths . 

2. A^A^ 4 X^E vertical deexcitation energy. 

The vertical transition energy was also calculated at the 
equilibrium geometry of the ^A-^ excited state (Table I), and the 
results are given in Tables V-VI . The Cl extrapolations used the same 
thresholds given in the previous section. In this case the inclusion 
of electron correlation causes an increase in the transition energy 
that is about twice as large as for vertical excitation. This is 
mostly because correlation of the C-0 sigma bonding electrons in the A E 
state becomes increasingly important at larger C-0 distances. This is 
reflected by the fact that in the E wave function the coefficient of 
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the ;j 4 tr * excitation 

. . . (3a 1 ) 2 (4a 1 ) 2 (6a 1 ) 2 (le) 4 (2e) 3 , (11) 

increases from about 0.06 at the ground state equilibrium geometry to 
about 0.09 at that of the excited state. The vibrational correction 
was calculated by including the zero point energies of ^ 1-3 for the A-^ 
state and of v ^_2 for the 2 E state (Table III). 

The calculated transition energy of 3.52 eV can be compared with 

/vO 

the intensity of emission bands originating in the A^A^ (v'= 0 ) state. 
The work of Inoue et al .^ 13 * 3 gives the maximum emission intensity at the 
v"=3 band with an energy of 3.53 eV. The results of Sutoh et al . 3 are 
also in excellent agreement with the present calculations. For the 
A^A^Cv'sO) -> X^E(v -# ) progression in emission the Franck-Condon factors 

were also calculated, and the relative per cent values obtained for 
v‘* = 0-8 levels are: 8 , 35, 74, 100, 97, 70, 38, 15, and 4. The 
experimental data for this emission (Fig. 7 of Ref. 10b) is in 
relatively good agreement with these values. 

"V v 

3. The A-X system origin. 

The energy of the (0-0) band of methoxy is obtained by comparing 
the calculations for the two states at their respective equilibrium 
geometries. The vibrational correction is obtained by including all 
three a^ modes ^^-3 of each state. The result, as given in Table VI, 
is a calculated system origin of 3.96 eV . This number is to be 
compared with the most recent experimental assignment of the origin in 
both absorption 4 and emission ^ 13 to the band at 3.91 eV. This level of 
agreement is consistent with the quality of the theoretical results, 
with the almost exact agreement noted in the previous section being 
somewhat fortuitous. There has been some difficulty in determining the 
system origin of methoxy, with earlier reports ^ 3 '^’^ 3 assigning it to 
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spectral features at 3.78 and 4.07 eV. It is expected that the present 
theoretical procedures should estimate the transition energies with an 
uncertainty of 0.05-0.1 eV, placing the assignment of 4.07 eV slightly 
outside of the error bounds and the one at 3.78 clearly beyond them. 

The present results confirm the most recent experimental assignment 
<3.91 eV) as being correct. 

4. Use of two-configuration E reference state. 

As noted above, the configuration that introduces left-right 
correlation of the C-0 cr-bonding electron pair CEq. (11)], becomes 
increasingly important in the X E wave function at larger C-0 distances. 
At dissociation this configuration should be included in an MCSCF 
treatment and in the Cl root set in order to obtain a qualitatively 
correct description of the wave function. Since the C-0 distance is 
quite large at the A state equilibrium geometry, the adequacy of the 
single configuration approach to the present problem was assessed by 
repeating the ^E calculations with Conf. (11) included in the MCSCF 
wave function and in the Cl root set. The ^A^ state does not require a 
comparable multiconfiguration treatment because Conf. (7) strongly 
dominates the wave function at the molecular conformations being 
considered. Because the C-0 bond in the A^ state is formally a 
one-electron bond, the single configuration treatment gives 
qualitatively correct dissociation to the neutral methyl and oxygen 
atom species. Additional root set configurations are only required at 
C-0 distances of at least 5 a Q in order to correctly describe the 
electronic state of the oxygen atom. 

The ^E calculations with Confs. (6) and (11) in both the MCSCF 
wave function and the Cl root set yielded extrapolated Cl energies 
approximately 3-6x10“^ h lower than did the single configuration 
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calculations. However, because some of the most important quadruple 
excitations, those that are doubly excited with respect to 
configuration (11), are explicitly included in the wave function, the 
Davidson correction aE^ becomes smaller in this multiconfiguration 
case. The net result of these two opposite effects is that the 
calculated system origin and vertical excitation energy would be 
decreased by 0.01 eV, and the vertical deexcitation energy would be 
increased by 0.04 eV. These changes would have no significant impact 
on the previous discussions. 

5. Analysis of the ^A^ wave function. 

A detailed analysis of both the E and ^A^ wave functions at the 
ground state equilibrium geometry was presented in paper I. The 
present work does not require any elaboration on that discussion but 

9 

makes possible an analysis of the A^ wave function at its equilibrium 
conformation. The wave function is very similar to the one previously 
discussed in I. Since the C-0 bond length is much longer at the ^A^ 
equilibrium geometry than at that of the ground state, one might expect 
that the charge distribution would begin to rearrange in a way that 
would resemble dissociation. In particular, the 5a^ molecular orbital, 
which is the C-0 cr-bonding orbital at short C-0 distances, becomes the 
non-bonding p-orbital on carbon at dissociation. However, at the A-^ 
equilibrium geometry the 5a^ orbital. Fig. 3, has not begun to localize 
on carbon but is still very much a bonding orbital with the largest 
part of its density on oxygen. Based upon the most complete A^ Cl 
wave function used in this study, there is a charge transfer of about 
0.5 electron from the hydrogen atoms to oxygen, with the carbon atom 
being essentially neutral. The electric dipole moment is calculated to 


be 3.8 Debye. 
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V. SURFACE CROSSINGS AND PREDISSOCIATION 

A. Experimental background. 

The well-studied hydroxyl radical has an electronic structure 

V 

similar to that o£ methoxy. The A state of OH shows evidence of strong 
predissociation at the v'=5 vibrational level, with weaker 
predissociation extending to the higher rotational levels of the v'=0 
state. D In the case of methoxy, the individual rotational lines 

have not been resolved. However, the experimental lif etimes^* 3 are 
consistent with the lower A-state vibrational levels (through v'=6) not 
being significantly predissociated . Emission has not been observed 
from the higher vibrational levels, allowing for the possibility that 
they may be short lived. The spacing of the vibrational levels in 
methoxy is much smaller than in hydroxy, and the nonradiative stability 
of the lower levels may simply arise from them being below the energy 
of the intersecting dissociative states.. On the other hand, it may be 
that the nonradiative mechanism is much weaker in the CH^O case. To 
help clarify this situation, we have obtained theoretical estimates of 

V o 

the energies at which the A A^ state intersects the potential surfaces 
of the 4 E, 4 A 2 * and ^^2 dissociative states. 

B. Theoretical method. 

Initially, attempts were made to define a limited, valence space 
Cl procedure for location of the surface intersections. However, 
balanced treatment of the four states was not obtained, and the 
selected Cl method with extrapolation was used instead. Preliminary 
calculations were carried out in the vicinity of the surface 
intersections to determine which configurations should be included in 

9 

the root sets. In the A^ case it is known that at dissociation, 
several root set configurations are needed to correctly describe the 
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state of the oxygen atom. Since the surface intersections occur at 
values of Rgg that are considerably larger than the equilibrium value, 
calculations were carried out to assess the importance of the secondary 
configurations. At the relevant conformations, Conf. <7) accounts for 
over SOU of the Cl wave function (based on c criterion), and the 

most important secondary configurations have coefficient magnitudes of 
about 0.05 or less. Preliminary calculations on the 4 E, 4 A 2 , ant * 2 A 2 
states indicated that secondary configurations were even less important 
in those cases. Therefore, throughout this phase of the study, single 
space configurations were used for the SCF wave functions and for the 
Cl root sets. 

The SCF orbitals for each state were obtained and were transformed 
into the corresponding ANO's. Cumulative selection with -c in the range 
20-50xl0"' 3 h was carried out on master configuration lists containing 
all single and double excitations. For the dissociative states the 
principal (root 3et) configurations are 

. . . (3a 1 ) 2 (4a 1 ) 2 (5a 1 ) 2 (6a 1 ) 1 (le) 4 (2e) 2 4 A 2 , 2 A 2 (12) 

. . . (3a 1 ) 2 (4a 1 ) 2 (5a 1 ) 1 (6a x ) 1 (le) 4 (2e> 3 4 E, (13) 

in which the 3a^ orbital is almost entirely °2s' 4a^ and le represent 
the three C-H <r-bonds, the 2e pair are nonbonding 0 2 p orbitals, 5a 1 is 
a cr-type orbital mostly (90-95*) on oxygen, and 6a^ is a or* orbital 
located mostly on carbon. When considered for the 2 A^ state, the main 
changes are that the 5a^ (double occupied) and 6a^ (unoccupied) 
orbitals are more equally delocalized over both carbon and oxygen. To 
illustrate these differences, the 5a^ and 6a^ MO's of the 4 E state and 
the 5a^ 110 of the A^ state are presented in Figs. 4-6 for a 
conformation near the 4 E - 2 A^ intersection. 

Near the potential surface intersections the state energies and 
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their spacings were found to be much more sensitive to small changes in 
Rqq than to changes in the other internal coordinates. Over the range 
of C-0 distances that include all three intersections, the optimum C-H 
distance for the A^ state changes by only 0.01 a Q , affecting the 
energy by about 0.1 mh. &CH wou ^^ change by about 0.03 a Q if optimized 
for the dissociative states, with an energy change of not more than 1.5 
mh. The situation is similar when optimization of the HC0 bond angle 
is considered. ®hc 0 wou ^^ vary by as much as 2° (.3 mh) for the ^A^ 
state and 5° <1.3 mh) for the repulsive states if full optimization 
were allowed. Because the energies are relatively insensitive to small 
changes in R^jj an< ^ ®HC0* was decided to optimize these coordinates 
at one C-0 distance for the A^ state and to use the resulting values 
unchanged for all of the states and conformations being studied. This 
reduces the apparent complexity of each potential surface intersection 
from a two dimensional surface to a single point. It is expected that 
this point will be reasonably close to the lowest point on the full- 
surface of intersection, because further optimization of Rqjj and S^q 
would lead to relatively small energy changes. With R co fixed at 
1.98 A, Rqjj and B^gg were optimized for the A^ state and found to be 
1.08 A and 97°. These values were used in all subsequent calculations. 

With Rgpj and 0^00 at the above values, the energies of the 

intersecting states were calculated as functions of Rgg • These 
energies, which included Cl extrapolation and the Davidson correction 
for higher excitations, were plotted, and the points of intersection 
were located graphically. When computing the intersection energy 

9 

relative to the E ground state, the theoretical results were shifted 
to accommodate the experimental value of the X-A system origin <3.91 
eV) since it is thought to be slightly more accurate than the present 
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theoretical value C3.96 eV) . 

C. Results and discussion. 

The intersection energies of the state with each disociative 

state are given in Table VII. It is seen that the lowest energy 
intersection is with the 4 A 2 s *-ate and is slightly above the v ' =7 
vibrational level (Table IV, ^A^, ^3 normal coord.). The other 
intersections (with 4 E and ^ 3 ) are close together and, by 
extrapolation of the data in Table IV, are located between the v'=9 and 
10 levels. Assuming an error in the intersection energy of no more 
than 0.1 eV, the corresponding vibrational quantum numbers are 
uncertain by no more than unity. These results suggest that the CH 3 O 
radical would not be strongly predissociated in the v'= 0-6 levels even 
if the nonradiative interaction were intrinsically strong. If the 
interaction is strong, the vibrational levels above v'=7 could be 
seriously affected, and emission from them might be quenched. In 
fact, there have been no experimental observations in emission of 
levels above v'= 6 . The absorption spectrum of Wendt and Hunziker 4 
contains a band identified as v'= 8 , but it has not been determined if 
individual rotational lines are broadened in either this or the 
lower-lying bands. 

It is interesting to compare the CH 3 O surface intersections with 
those of the hydroxyl radical. In the OH case, theory and 
experiment ' indicate that the the weak predissociation of the 

'VO 4. -v 

v'=0-2 levels of the A L state is due to the intersection of the A 
state potential curve with the lowest lying dissociative state( 4 E~). 

In the case of methoxy, the 4 A 2 3 tate is the lowest lying of the 
dissociative states and its electronic structure is analogous to that 
of the 4 IT state of hydroxy. A stronger predissociation occurs^® in 
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v oq v 

the hydroxy A state near the v'=5 level. Theory 0 indicates that the A 
state is crossed by both the 2 E - and 4 TT repulsive states in the 

9 — 

vicinity of the v'=5 level, with the E state being at slightly lower 

energy. The analogous CH 3 O dissociative states, 2 A 2 and 4 E, are also 

close together and cross the A state surface between v'=9 and 10. 

However, in this case, the ^2 state is predicted to lie slightly above 

the E state . If the comparison between the OH and CH 3 O radicals is 

valid, it would be reasonable to expect a weak predissociation in the 

CH 3 O ^Aji state near its intersection with the 4 A 2 state <v'=7), with a 

2 4 

stronger effect near the intersections with the A 2 and E states 
< v ' =9-10 ) . 

Simple considerations would lead one to predict that the 2 A^ - 4 E 
nonradiative interaction might be stronger than either the 2 A^ “ 4a 2 or 

9 9 

A^ - A 2 ones. Comparison of Conf. (7) with Confs. (12)-(13) shows 
that the principal 2 A^ and 4 E configurations differ formally by one 
spin orbital < 2 e -> 6 a^), whereas in both of the other cases the 
difference is two spin orbitals <2e 2 4 5a^6a^). The perturbation 
responsible for the doublet-quartet coupling interactions is the spin 
orbit operator, a sum of one-electron operators. The dominant 
contribution to the 2 A^ - 2 A 2 interaction is probably from the off 
diagonal nuclear kinetic energy operator, 

O?<q; 0 > I 3/3Q I (q;Q) > ( B/3Q) , 

where 'I' and are electronic wave functions for the two states, q 

represents the electronic coordinates, Q represents the nuclear 

coordinates, and the integration is over the electronic space. This 

4 1 

term can be shown to be equal to 

<¥<q:Q) I 3U/3Q1*' <q;Q)X3/3Q) 

E(Q)-E'CQ) * 
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where E and E' are the electronic state energies and 1A is the 
electrostatic potential energy operator. The dependence of U upon Q 
arises from the electron-nuclear attraction operator and, therefore, 
3U/9Q is a sum of one-electron operators. Thus, to the extent that the 
principal configurations C(7), (12), and (13)] dominate the full wave 

functions and that they are based on the same set of molecular 
orbitals, the rules 4 ^ for evaluating determinantal matrix elements 
would indicate that these (one-electron) perturbation operators would 
yield a nonzero coupling only for the ^A^ - 4 E case. 

In actuality, of course, the orbitals for the various states are 
not identical. Note in paticular the difference between the 5a^ 
orbitals in Figs. 4 and 6. Also, there are many important secondary 
terms in the full wave functions. Although the above discussion cannot 
be applied rigorously, it does make it plausible that the ^A^ - 4 E 
interaction would be considerably stronger than the other two. This 
would result in the strongest predissociation occurring near the A^ - 
4 E surface intersection. While this has not yet been observed in the 
methoxy case, the same general argument can also be applied to the OH 
radical and is consistent with the observations that have been made on 


that radical. 
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VI. SUMMARY AND CONCLUSIONS 
A. General. 

o 

The present study has resulted in a characterization of the 
excited state of CH 3 O that is consistent with the usual picture of it 
as arising through the excitation of an electron from the C-0 cr-bonding 
orbital into a nonbonding oxygen p-orbital. This results in the 
excited state having a much elongated C-0 bond, a relatively shallow 
potential well, and a large electric dipole moment. Accurate 

-v ^ 

transition energies have been calculated for the A-X system and agree 
very well with spectroscopic data. Harmonic frequencies for the three 
totally symmetric vibrational modes were calculated and were found to 
be in reasonable agreement with the available experimental assignments. 
By modeling the CH 3 O radical as a psuedo diatomic, it has been possible 
to obtain anharmonic vibrational levels for the C-0 stretch and to 
calculate Franck-Condon factors between the A and X states. These 
factors were found to be in reasonable agreement with the published 
experimental band intensities. 

In general, the present study agrees with and confirms the most 
recent spectroscopic assignments arising from both the absorption 4 and 
enisaion^^. Although in the past there has been some uncertainty in 
assigning the A-X system origin, this study confirms that the most 
recent experimental value of this quantity is correct. It is also 
predicted that predissociation of the methoxy radical into an oxygen 
atom and methyl radical would require excitation into the v'=7 or 8 
level of the A state, corresponding to a minimum photon energy of 4.47 
eV < *=2774 A). While this has not been confirmed experimentally, it is 
consistent with the observations. Bands corresponding to the levels 
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below v ' =7 are present, in absorption and emission spectra. These bands 
do not have significantly shortened lifetimes and, based upon these 
calculations, are not expected to show the effects of predissociation. 
Bands corresponding to levels with v'=7 and 8 are identified in 
absorption , but it has not been determined if the individual 
rotational lines are broadened. 

B. Importance of methoxy photolysis. 

A question addressed by this study was whether or not the 
photodissociation of methoxy, Eq. (5), is an important enough process 
to warrant inclusion in atmospheric models. First, considering Eq. (5) 
as a source of methyl radicals, it must be compared to other methyl 
sources usually included in the models. Three reactions generally 
considered to be important in this regard are: 


ch 4 ■ 

► OH 

-» ch 3 ♦ 

O 

N 

X 

(1) 

ch 4 ♦ 

Q 

\y 

o 

■* ch 3 

♦ OH 

(14) 

ch 4 

-*• Cl 

4 CH 3 ♦ 

HC1 

(15) 


Using recommended rate constants 43 and concentration profiles for the 
various species* - ^, the rates of methyl production due to the three 
reactions were estimated at several altitudes. Over the altitude range 
0-40 km these reactions always provide a rate of methyl production in 
excess of 10 molecules cm s . At most altitudes, Rxn . <1) is the 

most important CH 3 source, although Rxn. (15) is of comparable 
importance at 20 km and Rxn. (14) becomes significant at higher 
altitudes. The methyl production rates and data used in making these 
estimates are summarized in Table VIII for Rxn. (1). 

To estimate the rate of methoxy photolysis, several approximations 
have been made. Since we wish to know whether or not it is acceptable 
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to omit Rxn. (5) from the models, each approximation haa been made in 
such a way ao as to maximize the apparent photolysis rate. Thus, any 
recommendation to omit Rxn. (5) will be conservative in nature. 

In choosing a wavelength range for the 3 olar flux, the long 
wavelength limit is taken to correspond to the minimum energy (4.47 eV) 
necessary for predissociation as determined in this study. The short 
wavelength limit is chosen to correspond to the energy of the highest 
lying vibrational level of the A^ state into which excitation will 
occur. The diminishing magnitude of the Franck-Condon factors, as 
given in Section IV.B.l, supports the omission of levels above v'=12. 

By graphical extrapolation of the data in Table IV ( A^ , 'ig normal 
mode), the v'=12 level is estimated to be at 7274 cm -1 , requiring an 
excitation energy from the ground state of 4.77 eV. From the available 
tables of solar flux data 44 , we have chosen the wavelength range 
2778-2597 A as closely approximating the desired range of excitation 
energies, 4.47-4.77 eV. 

The absorption cross section 45 <r is taken to be 1.9 x 10 -1 ^ 
cm^/molecule over the entire wavelength range. This corresponds to the 
maximum value of the decadic extinction coefficient <50 1 mole -1 cm -1 ) 
found in the absorption spectrum 4 . The predissociation efficiency <t> is 
taken as unity over this wavelength range. These values of <r and <p are 
obvious over-estimates, but should provide an upper limit for the 
predissociation rate. The product J <r<p is an effective first order rate 
constant for CH 3 O predissociation and its calculation is outlined for 
several altitudes in Table IX. Using model 44 CH 3 O number densities, 
the calculated rates of photolysis are also presented in this table. 
Comparison of the photolysis rates in Table IX with the chemical rate 
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estimates in Table VIII shows that predissociation is unimportant as a 
methyl source throughout the troposphere and stratosphere. 

Consequently, its effect on the odd oxygen balance of the atmosphere is 
negligible at all altitudes. Even at the highest altitudes considered, 
unrealistically high CH 3 O concentrations of 10^ cm'^ would be necessary 
to provide a significant source of methyl radicals. 

The solar fluxes listed in Table IX vary strongly with altitude 
because the wavelength range used falls entirely within the absorption 
of the ozone layer. If the present calculations had indicated that 
predissociation could occur at longer wavelengths, the solar flux 
values used would be much greater. For example, if we were to include 
all methoxy vibrational levels v"= 0-12 in the predissociation 

range, the long wavelength limit would become equal to the system 
origin <3.91 eV or 3171 A), and the solar flux at all altitudes would 
be approximately 10^ photons ci’^ s - ^. Throughout the altitude range 
in Table IX, this flux would still produce a negligible rate of methyl 
production . 

At altitudes less than 10 km, the CH 3 O concentration may be 

considerably larger than the values given in Table IX. If CH 3 O is 

considered to be in a steady state due to Rxns. (3) and (4), its 

concentration is a function of the rate constants k 3 and k 4 and of the 

9 

03 , NO, and CH 3 O 2 concentrations. The model study of Logan et al. 

Q _Q 

suggests maximum CH 3 O 2 and NO concentrations of 5x10° molecules cm 
and approximately 100 ppt, respectively. Using these concentrations 
and current values of the rate constants 4 ^, a maximum CH 3 O 
concentration of 1.5x10^ molecules cm'^ is obtained. Combined with a 
maximum solar flux of 10 photons cm s A , 


this yields a photolysis 
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rate of roughly 0.3 molecules cm’^ s _1 . There are enough "order of 
magnitude" numbers in this rate that it would warrant further 
consideration if the flux could realistically be that large. However, 
the present calculations are accurate enough to make us confident that 
such low lying vibrational levels are not predissociated. 

Other photodissociation processes, such as direct excitation from 
the ground state of methoxy into its repulsive states or dissociation 
to the 0<^D) asymptote, have not been considered here. However, the 

.'W' 

semiquantitative potential curves of Fig. 1 suggest that the A - X 
vertical excitation is the only one at reasonably low energies. Direct 
dissociation to OC^D), requires an excitation wavelength less than 2130 
A and could not be important below 40 km. 

a/ •%/ 

One additional possibility is that the A f- X excitation of methoxy 
would be followed by interconversion to upper vibrational levels of the 
ground state, with subsequent dissociation. The thermodynamic CH 3 — 0 
dissociation energy 4 ^ <89 kcal/mole) and the calculated 19 barrier for 
H — CH 2 O dissociation <34 kcal/mole) are both less than the energy of 

/*V 'V 

the A - X system origin <3.91 eV, 90.2 kcal/mole). However, this 
interconversion is not expected to be fast because of the large A^ - 

9 

E "energy gap", a situation that is similar to the generally slow SI 4 
SO interconveraion of singlet states. Furthermore, as noted 
previously, the lifetimes of the A^ levels with v' less than 7 have 
been measured and do not show evidence of this dissociation 1 ®* 5 . 

Rxn, <5) also serves as a sink for the CH 3 O radical. However, the 
dominant sink is hydrogen abstraction by O 2 to form formaldehyde, Rxn. 

< 4 ) . Due to the abundance of 03 . the rate of this reaction is so great 
that photolysis is completely negligible as a CH 3 O sink <by at least 
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six orders of magnitude) . 
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TABLE I. Theoretical Values of the Ground and Excited State Equilibrium Geometries of the 
Methoxy Radical. 


Method 

R co 0(A) 

R CH e < A) 

0 HCO S(de9) 

Ref. 

2 E State 





First-Order Cl 

1.413(4) 

1.112(4) 

109.0(5) 

this study' 

MBPT 

1.41 

• • • 

• • • 

17c 

GVB-CI 

1.410 

1.111,1.112 

106.9,111.1 

18 

UHF 

1.386 

1.087 

109.9 

19 

2 A X State 





First-Order Cl 

1.614(6) 

1.091(2) 

102.6(2) 

this 3tudy' 

SCF 

1.65 

1.08 

102 

15 


a. Numbers in parenthesis are uncertainties in the last figure due to precision limitations 
of the potential surface analytical fit and of the search for its minimum. These values do 
not reflect the inherent accuracy limitations of the present Cl method. 
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TABLE II. Quadratic Force Constants for the Methoxy Radical. 



Force Constant® 

X 2 E State b 

'A 2 A 1 State b 


k 200 

k 110 

k 020 

k 101 

k 011 

k 002 


0.163 



0.968 

X 

10" 1 

0.454 

X 

10" 1 

0.715 

X 

10" 1 

0.527 



0.597 



0.278 

X 

10‘ 2 

0.333 

X 

10' 2 

0.226 

X 

10' 3 

-0.671 

X 

10~ 3 

0.122 

X 

10" 3 

0.115 

X 

10‘ 3 


a. The symbols k^ 3k are defined in Eq(9) of the text. 

b. The constants k 110' and k 020 are in hartree/bohr^; k 1Q1 and k Q11 are in 

O 

hartree/bohr-deg; and k Q Q2 ia in hartree/deg . 
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TABLE III. Vibrational Energies (in ca”*) for the CH3O and CD3O Radicals. 


Vibrational 

Mode 


CH 3 0 


CD 3 0 



This work® 

HF/631-G° 

Exptl 

This work® 

Exptl 

^E State 


3071 

3188 

. . . 

2203 

a a a 

* 2 

1335 

1585 

1325 c 

955 

1010 c 

*3 

1048 

1130 

1022 d 

1071 

1005 d 

^A^ State 

h 

3239 

• • • 

a a a 

2298 

a a a 

*2 

1308 

a a a 

a a a 

968 

a a a 


733 

a a a 

676 e 

732 

644 d 


a. calc ham frequencies. 

b. Ref. 19 

c. Ref. 13 

d. Ref. 10b 


e. Ref. 4 
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TABLE IV. Calculated Vibrational Energy Levels of Methoxy as a Pseudo Diatonic Molecule. 


State 

Vibrational 

Coordinate 

a - 1 

Energies 0 (cm ) 


2 E 

C-0 distance** 

544, 1623, 2685, 3729, 4753, 
6733, 7682, 8598, 9472, ... 

5755, 


normal mode c 

536, 1600, 2648, 3678, 4688, 
6640, 7577, 8481, 9249, ... 

5676, 

% 

C-0 distance** 

375, 1119, 1853, 2575, 3285, 
4661, 5326, 5973, 6602, ... 

3981, 


^3 normal mode c 

352, 1048, 1729, 2393, 3038, 
4256, 4821, 5361, 5884, ... 

3660, 


a. Relative to the minimum of the potential curve. 

b. Distance from oxygen atom to the center of mass of methyl moiety. 

c. As determined from three-dimensional harmonic analysis. 
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TABLE V. 

Calculated Energies 

(in hartrees) . 




State 

SCF 

CI(x=15mh) 

Cl ( x40) 



•V 


X state 

eq . conf oraation 




2 E 

-114.44560 

-114.72775 

-114.73849 

-0.02476 

K) 

*-* 

-114.29285 

-114.57127 

-114.58231 

-0.02386 

-a/ 

A state 

eq. conf oraation 




2 E 

-114.42331 

-114.70766 

-114.71852 

-0.02625 

2 », 

-114.30483 

-114.58179 

-114.59280 

-0.02445 



Page 42 


TABLE VI. Transition Energies (in eV) for the A 2 A^ - X 2 E System of the Methoxy Radical. 


Level of 
Theory 

Vert. Excitation 3 
2 A x f 2 E 

Vert. Deexcitation^ 

System Origin 0 
(0,0) band 

SCF 

4.16 

3.22 

3. S3 

Cl (t=15mh) 

4.26 

3.42 

3.97 

Extrap. Cl 

4.25 

3.42 

3.96 

Extrap. Cl ♦ aE^ 

4.27 

3.47 

3.97 

Extrap. Cl ♦ aE^ 
*■ vib corr 

4.21 

3.52 

3.96 


a. Both state energies calculated at ^E eq. conformation. 

b. Both state energies calculated at 2 Aj^ eq. conformation. 

c. Each state energy calculated at its own eq. conformation. 
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^9 

TABLE VII. Energies of Intersection Between the rAj and Repulsive States of the Methoxy 
Radical . 


Intersection C-0 distance 

at intersection 

(A) 


Energy relative 
to ^A^ minimum 
<c* _1 ) 


Energy relative 
to ground state 
(eV) 



2.01 

2.07 


4883 

6222 


4.47 

4.64 


2.09 


6398 


4.66 
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TABLE VIII. Rate of aethyl production due to reaction of the hydroxyl radical with 
nethane. 


Altitude 

(kn) 

•J*fl 

(K) 

rate const. ** k 
(ca 3 /aolecule s) 

tCH 4 ]C 3 

(aolecule/ca > 

COH] d 

(aolecule/ca 3 ) 

rate e 

(molecule/c® s) 

0 

291 

6.7xl0 -15 

3.5xl0 13 

10 6 

2.3xl0 5 

20 

219 

9.8xl0" 16 

1.7xl0 12 

10 6 

1.7xl0 3 

40 

268 

4.1xl0“ 15 

2.4xl0 10 

3xl0 7 

3.0xl0 3 


a. Ref. lc 

b. k=2.4xl0' 12 exp(-1710/T) for OH ♦ CH 4 4 CH 3 ♦ H 2 0, Ref. 43. 

c. derived from Fig. 8, Ref. la 

d. Fig 27, Ref. 2a and Fig 14a, Ref. la 

e. rate=kCCH 4 ] COH] 
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TABLE IX. Calculated photolysis rate of the methoxy radical. 


Altitude 

(km) 

Solar Flux 
J a 

2 

(photona/cm^ a) 

Photol . Const . 
k b 

CCH3O3 c 

Q 

(molecule/cm' > ) 

rate** 

(molecule/cm 3 s) 

10 

l.SxlO' 8 

2.8xl0" 27 

28 

8.0xl0' 26 

16 

6.6xl0" 7 

1.3xl0 -25 

7.6 

9.9xl0" 25 

22 

6.4X10" 1 

1.2xl0" 19 

117 

2. lxlO -18 

28 

2.7xl0 7 

5.1xl0" 12 

62 

3.2xl0“ 10 

34 

5.6X10 11 

l.lxlO' 7 

135 

l.SxlO' 5 

40 

5.7xl0 13 

l.lxlO' 5 

162 

1.8xl0" 3 

46 

2.7xl0 14 

S.lxlO" 5 

121 

6.2xl0~ 3 


a. From NASA model (Ref. 44), integrated over range a = 277.8-259.7 nm. 

.IQ 2 

b. Calculated as equal to J<r«, with <r=1.9xl0 c» / molecule and p=l molecule/photon. 

c. NASA model (Ref. 44), diurnal average, and at photochemical equilibrium. 

d. rate=k phot CCH3O] . 
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FIGURE CAPTIONS 


Fig. 1. 


Fig. 2. 


Fig . 3 . 


Fig . 4 . 


Fig . 5 . 


Fig . 6 . 


Potential energy curves of the methoxy radical at C 3v 

conformations. R,-.^ = 2.041 a Q ; B^CO decreases from 113 to 

90° as R^q increases from 2.27 to 5.67 a Q . (Reproduced from 

Ref. 3) 
o 

A^ Potential energy curves of the methoxy radical as 
functions of the C-0 distance and ^3 normal coordinate. aR C0 
is obtained in bohra by multiplying the abscissa value by 
0.378 (dashed curve) or 0.359 (solid curve). Energy levels 
are for the dashed curve. 

The 5a^ SCF molecular orbital for the CH 3 O 2 A^ state, in an 
HCO plane, at the 2 A 1 equilibrium conformation: R co = 1.61 A, 
Rch = 1*0^ A, and 9 hC 0 = 102.6°. Contours enclose 90, 70, 50, 
30, and 10* of integrated three dimensional probability. 

Major tic spacing = 2 a Q . 

The 5a^ SCF molecular orbital for the CH 3 O ^E state, in an 
HCO plane, at the conformation: R^q = 2.08 A, Rq^ = 1.08 A, 

and BhC 0 = 97* 0 °- Contours and tic spacing as in Fig. 3. 

The 6 a^ SCF molecular orbital for the CH 3 O 4 E state. 
Conformation, contours, and tic spacing as in Fig. 4. 

The 5a^ SCF molecular orbital for the CH 3 O 2 A^ state. 
Conformation, contours, and tic spacing as in Fig. 4. 




to 

QJ 

5 


03 


5 s 

<b 

>2 


-114.00 


-114.10 


-114.20 


-114.30 


-114.40 


-114.50 














